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Abstract 

We generalize the Arbitrage Pricing Theory (APT) to include the contri- 
bution of virtual arbitrage opportunities. We model the arbitrage return by a 
stochastic process. The latter is incorporated in the APT framework to calcu- 
late the correction to the APT due to the virtual arbitrage opportunities. The 
resulting relations reduce to the APT for an infinitely fast market reaction or in 
the case where the virtual arbitrage is absent. Corrections to the Capital Asset 
Pricing Model (CAPM) are also derived. 

1 Introduction 

Along with the Capital Asset Pricing Model (CAPM) [0, |^, ^ the Arbitrage Pricing 
Theory (APT) 0] is a standard model for determining the expected rate of return on 
individual stocks and on a portfolio of stocks |^. In contrast to the CAPM, the APT 
does not rely on any assumptions about utility functions or the assumption that agents 
consider only the mean and variance of possible portfolios. Taking into account criti- 
cal remarks about the expected utility approach to the theory of the decision making 
under uncertainty, and the long running discussion about the definition of risk, this 
seems to be a considerable achievement. What the APT does imply are homogeneous 
expectations, assumptions of linear factor weightings, a large enough number of secu- 
rities to eliminate the specific risk and no-arbitrage condition. In simplified terms, the 
latter means that the return on a riskless portfolio should be equal to the riskless rate 
of return which can be taken for a moment equal to the rate of the return on a bank 
deposit (here, for sake of simplicity, we assume perfect capital market conditions, i.e. 
single rate of lending and borrowing without restrictions, absence of transaction costs 
and bid-ask spread; we touch on possible generalization in conclusion). Being a basis 
for the APT and derivative pricing P], the no-arbitrage assumption appears to be very 
reasonable and robust. Indeed, if any arbitrage possibility would exist then agents (ar- 
bitrageurs) would use the opportunity to make an abnormal riskless profit which itself 
will bring the system to the equilibrium and eliminate the arbitrage opportunity. Thus 
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the arbitrage cannot exist for long and its lifetime depends on a liquidity of the market. 
This, however, does not mean that the arbitrage opportunities do not exist at all and 
cannot influence the asset pricing violating the APT assumption. Moreover, there exist 
empirical studies which point out the existence of virtual arbitrage opportunities even 
for very liquid markets (order of 5 minutes for futures on S&P and can be much 
longer for less liquid markets such as bonds market and derivative market [Q). That is 
why in this paper we try to overcome the no-arbitrage assumption and suggest a model 
to account for the existence of virtual arbitrage opportunities and their influence on 
asset pricing in the framework of the APT. To this end we follow the line introduced 
in Ref. ^ to account the virtual arbitrage in derivative pricing. 

One of the possible approaches to the problem of handling of virtual arbitrage has 
been suggested in Ref [|l3 and based on a field-theoretical description of arbitrage 
opportunities and the corresponding money flows. Being a consistent theory, the ap- 
proach is however extremely complicated and does not allow one yet to get simple 
analytical results which would be easy understandable and handleable. That is why in 
this paper we develop simplified tractable analytical version to account for virtual ar- 
bitrage opportunities which does not use complicated techniques and results in simple 
enough final formulas. 

The paper is organized as follows. In next section we derive an effective equation 
for a price of a riskless portfolio in presence of virtual arbitrage opportunities. There 
it is shown how the local arbitrage opportunities can be introduced in the model and 
how they change the rate of return. In section 3 we derive equations of the APT under 
virtual arbitrage. These equations generalizes the APT relations and converge to them 
in the limit of absence of the arbitrage opportunities or infinitely fast market reaction. 
Section 4 is devoted to the corrections to the CAPM. In conclusion we discuss the 
drawbacks of the model and possible ways to improve it. 

2 Effective equation for riskless portfolio in pres- 
ence of virtual arbitrage 

In this section we introduce a model for virtual arbitrage fluctuations and describe 
their influence on return on riskless portfolio. Under no- arbitrage assumption the rate 
of return on a riskless portfolio shall be equal to the single riskless rate of return, which 
is equal to, for example, the return on a bank deposit. However, if we allow existence 
of virtual arbitrage opportunities, the rate of return on a portfolio does not have to be 
equal the rate of return on bank deposit, it is random (since unpredictable) and might 
depend on the structure of a portfolio. Indeed, if, for example, the portfolio consists of 
only riskless assets (bank deposit or bond) there should be no arbitrage opportunities 
with this portfolio and the rate of return on this portfolio, by definition, shall be equal 
to the rate of return on a bank deposit tq. In the general case of a riskless portfolio 
which consists of diversified risky assets the existence of the arbitrage fluctuations 
would depend on the portfolio composition. To give an intuitive example, we can 
imagine that if some market sector is more attractive for investors, the trades and 
resulting arbitrage fluctuations will occur more often. This will result in more arbitrage 



opportunities for portfolios containing the assets of the sector and the corresponding 
composition dependence. 

Let us consider riskless portfoho 11 which is created by A^ + 1 assets with fractions 
{xi}iLo- In the case of no-arbitrage the portfoho IT would satisfy the following equation: 

^-ron = o, n(o) = i (1) 

where tq is riskless interest rate on bank deposit. However, in case of the virtual 
arbitrage RHS of Eq(|l]) shall be changed to 

7^(t,^)^ 

where TZ{t, U) represents the virtual arbitrage return. To find an expression for TZ{t, U) 
let us imagine that at some moment of time r < t a fluctuation of the return (an 
arbitrage opportunity) appeared in the market. We then denote this instantaneous 
arbitrage return as z/(r, 11). Arbitragers would react to this circumstance and act 
in such a way that the arbitrage gradually disappears and the market returns to its 
equilibrium state, i.e. the absence of the arbitrage. For small enough fluctuations it is 
natural to assume that the arbitrage return TZ (in absence of other fluctuations) evolves 
according to the following equation: 

— = -A7^, 7^(r)=z/(r,^) (2) 

with some parameter A which is characteristic for the market. This parameter can be 
either estimated from a microscopic theory like |1^ or can be found from the market 



using an analogue of the fluctuation-dissipation theorem [jTT|. In the last case the 
parameter A can be estimated from the market data as 



(3) 
and may well be a function of time, portfolio composition and even prices of assets. 
In what follows we however consider A as a constant to get simple analytical formulas 
for APT corrections. The generalization to the case of time-dependent parameters is 
straightforward. 

The solution of Eqn(|^) gives us TZ{t, 11) = z/(r, n)e~'^^*~^'' which, after summing over 
all possible fluctuations, leads us to the following expression for the arbitrage return 

n{t,U)= f e^^(*-"V(r, n)dr . (4) 

J —OQ 

To specify the stochastic process i^(t,n) we assume that the fluctuations at different 
times are independent and form the white noise with a variance S^(n) which depends 
on the structure of the portfolio 11: 

(Kt,n)), = , {u{t,uMt',u)), = j:\u)-6{t-t') , (5) 



where, again, the parameter S^(n) can be taken from the market: 

S^(n)/2A = ((rn-ro)n__ . (6) 



market 

To simphfy the consideration we assume here that the quantity does not depend on 
time but this hmitation may be straightforwardly overcome (see discussion in the con- 
clusion) . 

Since we introduced the stochastic arbitrage return 7l{t, 11), Eqn(P has to be sub- 
stituted with the following equation: 

^-ro^ = 7^(t,^)^ (7) 

or, in the integral form, 

U= f G{t,t')n{t',U)U{t')dt' (8) 

where G{t, t') = 0(t — t')e^'"*^*~* ^ is Green function of the problem: 

(^ - ro)G(t,t') = 5it - t') , G(t,t')l*<t' = . 
We can iterate Eq(|^) and substitute Eq(^ in RHS of Eq(0) which gives 

— - roH = 7^(t, n) / G{t, t')n{t', U)U{t')dt' . (9) 

(It J —oo 

Our next step is to average Eq (H) over virtual arbitrage fluctuations and to obtain 
the effective pricing equation for the average price n. At the first order in 1/A it can 
be written as 

^-ron=/ G{t,t')K{t,t',n)n{t')dt' . (10) 

dt J —oo 

where the kernel 

ir(t,t^^) = (7^(t,^)7^(t',^)), 

is given by the expression 0: 

K{t, t') = ^^0{t - t')e-^(*-*') + ^^0it' - t)e-^(*'-*) (11) 

which can be easily obtained from Eq(^). Collecting everything together, Eqs(|TU|,|Tip 
are used in place of Eq(|l]) in situations where virtual arbitrage opportunities exist. 
To solve Eq(PIlD we first notice that its RHS can be written as 

/* Git,t')K{t,t',n)n{t')dt'= f e-^(*-*')e^°(*-*')^^II(t')(it' • 

J —oo J —CO Za 

We now make an approximation T[{t') = e^°* because the correction to this expression 
is order of ^^ ' and is irrelevant for our consideration of RHS of Eq(|lOD since it was 



derived within this order of accuracy. Such approximation results in the following 
relation: 

/ G{t, t')Kit, t', n)nit')dt' = ^r^ e^o* . 

If we substitute it into Eq (|iy) this gives us the following approximate differential 
equation for the average portfolio price: 



which has the solution: 



2A2 



n(t) = (1 + -^^t)e^ 
or, to the same level of accuracy, 

n(t) = e(''°+^)* . (13) 

The latter expression we use in the next section in place of the no-arbitrage expression 
n(t) = e^°* in the Arbitrage Pricing Theory. 

Let us emphasize that we, as well as an investor, are interested in the average value 
of the portfolio and the corresponding rate of return on the average portfolio rather 
than in average rate of return on the riskless portfolio which is exactly equal to tq 
and, in fact, is not influenced by the arbitrage. Indeed, the only valuable and material 
thing for an investor is the average value of her investment. She is not interested in a 
mathematical quantity which is not actually connected with her wealth. The difference 
between the two returns is the first effect which is solely due to the presence of the 
virtual arbitrage. 

3 Corrections to APT 

Let us first remind ourselves of the standard derivation of Arbitrage Pricing Theory 
relations. The first assumption below APT is the existence of M fluctuating leading 
parameters, or factors, {C,j}fLi which define the prices of all A^ + 1 assets such that 

and 

ri=fi + Y,bij^j + ^i (14) 

j 

with bij = {ri — fi,^j)^ as a measure of influence of j-th parameter on i-th asset and 
ej is a residual risk which is independent on the factors and can be eliminated in large 
portfolios. To form a riskless portfolio we pick up such fractions {x,} that 

Y^ Xibij = 



for any j = 1, ..., M and they obey normalization condition J^i^i = ^- We assume that 
0-th asset is riskless, e.g. bank deposit or treasury bonds. 

The return on the portfoho is J2i ^i^i = J2i ^i^i ^^^ has to be equal to the risk-free 
interest rate tq in the case of no-arbitrage. This leads to the equation 

J2riXi = ro (15) 

i 

or, after a change of variables Xi = yi + 6oi, 

i 

subject to the constraints 

J2 Vi^ij = V j , J2yi = . 

i i 

The solution of the equation for fj then can be written as 

fi = a + J^bijlj ■ (16) 

3 

where a and {7^} are set of independent parameters. Since the 0-th asset was chosen as 
risk-free, the coefficients boj = for all j and a = tq. Eq(|T^) constitutes the Arbitrage 
Pricing Theory relations. 

To generalize the APT relations to the case of virtual arbitrage we have to deal 
with average portfolios considered in the previous section. We define the rates of the 
return {r^}^o '^^ ^"^^ component of an average riskless portfolio as 

n(t) = eS«">'* . 

These are the quantities which have to substitute the rates of return in the APT 
relations. 

To find an expression for {f^}^Q we have to substitute tq in Eq(P^ on the return 
to the riskless portfolio derived under virtual arbitrage consideration in the previous 
section, i.e. on (rg H — 2a2^) ^^ follows from Eq(P^. This gives us the equation for f'f 

J:r'a' = ^ (17) 

subject to constraints 

El/A, = vj, E?/« = o- (18) 

i i 

To simplify the following consideration we introduce a convenient basis in the portfolio 
space such that the above constraints will take the form 

r?i = , < i < M 



in this new basis while {'>]j}f=M+i '^i^ be coordinates in the subspace of risk- free 

I A/ . 
ti=0- 

e' = 1 Vz , el, = 6„- Vz , 1 < J < M . 



portfohos. To this end we first introduce vectors {e'f}f^ 



These M+1 vectors are hnear independent (otherwise it would be linear dependence for 
the factors and M + 1 constraints (18) would not be independent). The rest of iV — M 



linear independent vectors {e^}^^.,.]^ can be chosen arbitrary but linear independent 
with {e^lf^Q. For example, if 

bij 7^ , Vj and < i < M + 1 

then one possible choice could be 

(i'i,j = ^ij yj ,M + i<t<N + i . 

The next step is the Gram-Schmidt orthogonalisation and normalisation of the vector 
set {e^}^o which gives us the basis set {eij^A/f+i- We shall start the orthogonalisation 
with the vector e'^: 





eo 


eo 




V/K>e^) ^ 


proceed with the vector e'^ as 






ei = ^ 




e'l - (e'i,eo)eo 




(e'l - (e'l, eo)eo, ei - (ei, eo)eo) 
and carry on following the standard procedure: 

^i = i^'j- 12i^P ^k)ek I (e;., e;.) - ^(e^, e^) 

\ fc=0 / \ fc=0 

It is easy to check now that in this new basis {ej}^Q the constraints ( |TB|) can be 
represented as a zero projection along the first M+1 basis vectors, i.e. as 

r^. = , < 2 < M 

where r^j are coordinates in the new basis. Indeed, for any vector y the conditions 

(y,e;) = 0, 0<j<M 

are exactly equivalent to the constraints (|l^). This remains true also after the orthog- 
onalisation 

(y,e,) = 0, 0<j<M 

due to step-by-step nature of the procedure. Therefore any vector 

Af+l 
j=Af+l 



represents a risk- free portfolio and the set {Gj}f^^.f_^i is a basis in the subspace of 
riskless portfohos. 

To express the riskless portfolios basis in terms of the original assets we introduce 
the rotation matrix U = \\Uik\\'- 



U 



ei 



(19) 



Using this matrix the portfolio represented by the basis vector e^ can be constructed 
by collecting j-th assets (j' = 0, .., iV + 1) with the fractions etj + 6oj. 

Now we are ready to discuss the possible form of S^(n) as a function of portfolio 
structure, i.e. as a function of {^ij^o ^^ {Vi}iLM+i- First of all we will not allow 
the existence of arbitrage fluctuations for a portfolio which consists of pure bonds, i.e. 
when all r]i equal to zero. This provides that the rate of return on the bank deposit is 
actually equal to tq, does not have any arbitrage corrections and our use of rg is self- 
consistent. It means that series expansion of S^(?/) cannot contain a constant term. 
Second, since the virtual arbitrage return is brought in by portfolio r], the arbitrage 
return on the portfolio —rj has to be equal to minus the latter and ^"^{f]) = E^(— 77). 
This results in the absence of all odd terms. Therefore the series expansion of S^(n) 
has the form: 

N 
iJ=M+l 

or 

N 

E\n)= J2 ^lVrVj + - , S2.fc = for 0<t,...,k<M. 

i,j=M+l 

In what follows we keep only first nontrivial term, [rj, S^r/), S^ = | |S^^| | in the expansion 
though the final result will be valid for the general case too. 

Returning back to pricing equation ( |T7|) we substitute y = Ur], S^(II) = (77,^^77), 
J2i fiiji = (f , y) we come to the pricing equation for r: 

(r, Ut]) = {t], Y?vi) 

subject to constraints 

r/i = , < J < M . 

The solution of the problem can be easily found as 

N ^2 

n = a + Y: h^i, + E U,,-f^r,i . (20) 

j i,k,l=0 ^^ 

At this point we face an unpleasant problem: a return on an asset in the portfolio 
depends on the structure of the portfolio. This is the second effect which stems from 
the virtual arbitrage presence. However unpleasant it is, it is hardly surprising since 
it was derived under the assumption of virtual arbitrage which is portfolio-dependent. 
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Intuitively it is also clear: if an asset constitutes a part of a "hot" arbitrage portfolio 
its rate of return will be different from assets in "quiet" portfolios ceteris paribus. 

To find an average growth factor for the z-th asset we have to take a sum of e'''-'^^' 
weighted with probabilities of appearance of the portfolio 11, i.e. to calculate the 
average value 

which defines the rate of growth fi of the average growth factor. This fj is a counterpart 
of the average rate of return under the virtual arbitrage presence since it characterizes 
how fast i-th asset grows in average. It is easy to see that in the first nontrivial order 

\ k,m,l,l'=0 

and, hence, with the same measure of accuracy the rate of growth is given by 

N v^2 ^2 

f, = a + J2k,^,+ 5: u,,U,„,^^{viVi')u (21) 

j k,m,l,l'=0 °^ 

Since the 0-th asset was chosen as risk-free, the coefficients Poj = for all j and 

j k,m,l,l'=0 °'^ k,m,l,l'=0 °'^ 

(22) 
The last equation looks quite complicated. It can be presented in a more compact form 
using matrix notations and the following matrix A: 

4A2 \ dr]i dr]k /„ ' 
The pricing relations then can be rewritten as 

2A2^ ' 2A 



ro + T. hn^ + —2^U-A- U-% - (U-A. f/-i)oo . (23) 



J 



Eqns (|22|j23|) gives a generalization of the APT relations in the case where virtual 
arbitrage opportunities exist. It is easy to see that the corrections are proportional to 
the square of the product of the variance of the arbitrage fluctuations and the square 
of a characteristic lifetime of the arbitrage. This results in the disappearance of the 
correlations in the case of infinitely fast market reaction (A — *■ cxd) or an absence of the 
arbitrage (S — > 0) and the reduction of Eqns( P^p5D to the APT expression (^B|). 



4 Simplest application: correction to CAPM 

In this short section we concentrate on a particular case of APT with the only one 
leading parameter which will be chosen to be a random part of a return on a market 
portfolio. In this case the standard APT relations are reduced to the CAPM equations. 

9 



In the case of only one leading parameter chosen to be a random part of a return 
on the market portfolio ^ = (r^ — fm)/am, Eqn ( p3D can be rewritten as 

fi=ro + hl^iU ■ A ■ U-% - ^{U ■ A ■ U-\o (24) 

where bi = {ri — fi^Tm — Trr) jom- The next step is to introduce the market portfolio 
fractions Qi and find an expression for the average rate of return on the market portfolio 
Ei^^ifi using Eqn (^: 

fm = To + Y.^Al + ^M™ - ^{U ■ A ■ f/-l)oo 

where we introduce a notation m™ = ^i{U ■ A ■ U^^)ii9i. This allow us to define a value 
of the variable 7: 

and to find a final expression for the CAPM generalized to the virtual arbitrage as- 
sumption: 

^^(- \^(^^ ((TT A rr-1^ ,„A , (^ ■ A ■ U'^h {U ■ A ■ U'^p 
ri = ro + f3i{r^-ro) + ^[{U-A-U )oo - u j + — — 

(25) 
with the standard definition Pi = {ri — r^, r^ — rm)/cr^. The first two terms represent 
the CAPM under the no-arbitrage assumption and the remaining terms are the virtual 
arbitrage corrections. Everything said in the last paragraph of the previous section 
can be repeated with respect to Eqn(P^ for the virtual arbitrage generalization of the 
CAPM. 



5 Discussion 

In conclusion, we introduced a virtual arbitrage opportunities in the framework of 
the Arbitrage Pricing and derived corrections to the APT pricing relations. These 
corrections disappear as E ^ or A -^ i.e. when there is no arbitrage fluctuations 
or the speed of market reaction to the mispricing is infinite (which corresponds to 
extremaly liquid market). 

In the course of the analysis we faced two major effects which are specific for the 
presence of the arbitrage. The first one is a necessity to analyze growth factors rather 
than the the corresponding rates of return. The second effect is a dependence of a 
growth factor of an asset on the structure of a riskless portfolio which includes the 
asset. This forced us to introduce an average over all riskless portfolios and reduced 
the problem to a calculation of the matrix A 
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which is essentially a new element. We already mentioned that S^(n) can be obtained 
from the market using Eqn(^. It means that A, in principle, can also be obtained from 
the market following (^61) . However the most economical procedure of the evaluation 
still has to be worked out. It might, for example, be reasonable to consider some most 
frequent arbitrages and reduce the space of the relevant riskless portfolios. Another 
option is to study thhe main factors generating arbitrage opportunities and carry out 
a factor analysis similar to the original APT. At the moment we leave this as an 
open question which requires further investigation. Now let us discuss some obvious 
drawbacks of the model and ways to improve it. 

First of all, we have to admit immediately that it would be more difficult to carry 
out empirical tests of the Virtual Arbitrage Pricing Theory than the APT because 
the additional terms in Eqn(0) do not make life easier in any respect while even 
empirics for the CAPM and the APT [|T^ have not produced yet any decisive result. 
Furthermore, all critical comments of APT analysis can be forwarded to the presented 
model, except for the no- arbitrage constraint. Homogeneous expectations, absence 
of clear recipe of how to measure the factors ^j and possible time dependence of the 
parameters of the model are three obvious points to criticize. These are not new 
problems and many efforts to overcome these difficulties have been undertaken. It is 
possible to demonstrate that the virtual arbitrage model can be improved in the same 
manner by these methods as they succeed for no-arbitrage APT analysis. This, in 
principle, allows one to include the transaction costs, taxes and market imperfection in 
the present model by rederiving the "bare" APT relations and then to add the virtual 
arbitrage in the consideration with minor modifications. 

The second point concerns the market reaction to the arbitrage opportunity, or, 
qualitatively the form of TZ{t, H) in Eqn(^. It may be argued that the market reaction 
is not exponential as assumed in Eqn(|^), but has another functional dependence. This 
dependence can be found from statistical analysis of the stock prices and then included 
in the equation for TZ{t, H) (in particular, the functional dependence can change with 
time, for example A in (0) can be a function of r). It certainly complicates the model 
but leaves the general framework intact. 

Another point to consider is the absence of correlations between virtual arbitrage 
opportunities which we assumed in the text, i.e. the white noise character of the process 
z/(t,H). It is clear that some correlations can be easily included in the model by a 
proper substitution of the relations in Eqn(^). Such generalizations, though making 
the analytical study almost impossible, allow one to proceed with numerical analysis 
for the model. 

Finally, the model contains new parameters such as E and A defined by Eqns (|^) and 
(0) from the market data. It might appear that time-independent parameters are not 
good enough approximation and both S and A are functions of time as a consequence 
of market intermittent busts of activity. Furthermore, the parameter A might actually 
depends on a structure of the portfolio exactly in the same manner as S does because 
of particular sector preferences of arbitrageurs. Both situations can be processed in 
a similar manner to the simplest case we considered above but it will make the final 
results much more complicated. 
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